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STRUCTURAL DYNAMICS OF THE HELIOGYRO 
SUMMARY 

The Heliogyro is a particular form of the solar sail in which the reflecting 
surface consists of a set of long, narrow and very thin blades that rotate about 
a central axis in the manner of a helicopter rotor. Blade pitch is the primary 
f means of control. The structural dynamics of the Heliogyro are studied in this 

i 

report in order to provide a basis for engineering development and to answer 
questions regarding concept feasibility. The topics that are examined include 
vibration modes, damping, coupling between blade pitch and blade bending, and 
dynamic response to maneuver commands. The results of a brief experimental 
program are presented and compared with analysis. 


INTRODUCTION 

The Heliogyro is a particular form of the solar sail, which is the name 
that has been given (Refs. 1, 2, 3) to a class of devices which develop useful 
mechanical force in space from solar radiation pressure. The principle of 
operation of the solar sail is illustrated in Figure 1. Light is partly 
reflected and partly absorbed at the surface of an object. The part that is 
reflected exerts a pressure 


p = p f cos 2 y (l) 

r K o r 1 

normal to the surface, where p is the total pressure for normal incidence due 

o r 

to reflected light, f is the coefficient of reflectivity, and y is the incidence 
angle. The part that is absorbed exerts a pressure 

p a = T P o (1 " V COS Y (2) 

in the direction of the incident radiation. The part that is absorbed is 
eventually reradiated diffusely from both sides of the reflecting surface. 


The components of radiation pressure are combined as 
pressure, p^, in the direction of illumination and a 
dicular to the illumination. 


follows to obtain a drag 
lift pressure, p^, perpen- 


= P 0 (f r cos 3 y + j ( 1 ”f p) cos y) 


(3) 


p £ = p Q (f r cos 2 y sin y) 


W 


The coefficient of reflectivity, f , is high for the metalized surfaces 
envisaged in the solar sail concept. It is commonly assumed equal to unity in 
preliminary design studies. 

The lift and drag components of solar radiation pressure have a strong 
analogy with the lift and drag components of aerodynamic pressure, so that a 
pseudo-aerodynamic theory of solar sail behaviour can be developed. The lift 
force is particularly important because it provides the means for generating 
forces to accelerate or decelerate a body that is in a circular orbit around the 
sun . 


The magnitude of the solar radiation pressure is not large by terrestrial 
standards. At the Earth's distance from the sun its average is equal to 
0.90 x 10 ^ dynes/cm 2 = 0. 188 x 10 ^ lb/ft 2 = . 1 304 x 10 ^ psi. Very large 
surface areas are required to obtain enough thrust to be useful for propulsion. 
Fortunately for the solar sail concept, the technology is available for the 
manufacture of extremely thin reflecting surfaces. Standard commercially avail- 
able l/A mil aluminized Mylar, for example, makes a solar sail that is quite 
competetive with other means of propulsion on long missions. 


The fundamental performance parameter of the solar sail is the sail light- 
ness number, X , which is defined as the ratio of the solar radiation pressure 
s 

to the solar gravitational force per unit area. Since both factors vary in- 
versely as the square of the distance from the sun, the lightness number is con- 
stant within the solar system. The sail lightness number is related to the 
thickness, h, of the sail, and its weight density, pg, by 


- 2.15 x 10 
s (pg) h 


( 5 ) 


2 


where the units of (pg)h are lbs/in 2 . The value of for 1/1* mil aluminized 
Mylar is approximately equal to 0.16. 

Figure 2, taken from Ref. A, shows a comparison of the propulsive efficiency 
of the solar sail and several other propulsion systems. Specific impulse is 
used as the measure of efficiency. It is seen that the specific impulse of the • 
solar sail increases (linearly) with mission duration and that it outstrips all 
competing systems for mission durations greater than 1000 days. The range of 
values of lightness number used in Figure 2, .25 < < 1.0, extends from a 

value that is available with present technology to a value that seems feasible 
with future development. The performance ranges for the other propulsion system 
were selected in the same spirit. 

Analyses of optimum interplanetary flight trajectories have been reported 
by several solar sail investigators (Refs. 5, 6, 7, 8), particularly for voyages 
to Mars and Jupiter. Applications of the solar sail concept for other purposes, 
including orientation control, station keeping, station changing, and micro- 
meteoroid collection, have also been studied, Ref. 6 and 9. 

In spite of favorable performance estimates, the solar sail was not 
initially given serious consideration, due to the absence of credible solutions 
for deploying, rigidizing and orienting the extremely large surface areas that 
are required. Table I shows a morphology of some of the choices that are avail- 
able in the design of the solar sail. Early investigators favored the circular 
shape, which virtually defies deployment in large sizes. The more recent Helio- 
gyro concept (Ref. 9), which employs the design choices underlined in Table I 
( very slender , centrifugal force , reflecting surface rotation) shows promise of 
providing practical solutions for the problems associated with large size. 

The Heliogyro is, in effect, a helicopter rotor operating in a space 
environment. Its blades are made from thin reflecting sheets that are stowed on 
spools prior to deployment. The rotation of reflecting surfaces is achieved by 
pitching the blades about their lengthwise axes. Centrifugal force provides 
effective bending stiffness and torsional stiffness for the blades. A sketch 
of a two-bladed experimental model, proposed in Ref. 9, is shown in Figure 3- 
Concept models of much larger multi-bladed vehicles for manned interplanetary 
voyages are also described in Ref. 9. 
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The design concept of the Heliogyro borrows much from helicopter technology 
including the use of centrifugal force to rigidize the blades, and the use of 
blade' pitch for control. The use of cyclic and collective pitch controls to 
produce lateral forces and moments and to impart initial spin during deployment 
are discussed in Ref. 9. Calculations are made there of deployment times, 
precession rates, and of the times required to perform such useful maneuvers 
as planetary escape. The results show that the maneuver performance is good 
and that the stress levels are low for vehicles ranging in size from the smallest 
up to several hundred thousand pounds gross weight, even though the latter may 
have rotor diameters of fifty miles or more. The calculations are, however, 
based on the usual assumption made in helicopter performance analysis that 
structural dynamic effects are negligible. Thus the rotor blades are assumed 
to respond instantly to pitch inputs, and pitching moments due to flapwise 
and chordwise motions are ignored even though the effective torsional stiffness 
provided by centrifugal force is very low. The means for providing damping in 
the absence of an aerodynamic environment were not discussed. 

The purpose of the present report is to treat the structural dynamics of 
the Heliogyro in such a manner as to provide a basis for engineering develop- 
ment and to answer questions regarding concept feasibility that are raised by 
dynamic considerations. The subject will be treated by the techniques of heli- 
copter rotor dynamics, but differences between an Earthbound and a space envi ron- 
ment, and differences in scale and proportion require that the treatment start 
at an elementary level and that it carefully examine the basis of every assump- 
tion. I t wi 1 1 be convenient, for the estimation of the magnitude of effects, 
to refer occasionally to the numerical values of parameters for a specific 
representative blade design. That design, whose properties are listed in Table 
2, will be called the Sample Blade. 
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UNCOUPLED EQUATIONS OF MOTION 


Since helicopter rotor blades are relatively slender (a condition which 
applies, a fortiori, to the Heliogyro blade), it is usual to analyze their 
motion by means of engineering beam theory. The notational system that will 
be used to describe blade motions is shown in Figure l t. The equilibrium 
position of the blade is assumed to lie in a plane perpendicular to the axis 
of rotation. The undeflected axis of the blade is assumed to be straight. 

The assumption of initial straightness, which is made for convenience in 
analysis, can be relaxed without difficulty. The independent degrees of 
freedom of the blade at points along its axis are its translations (u, v, and 
w) , and its rotation, 0, about the blade axis. Distortion of the cross- 
sectional shape (e.g., by cambering) is neglected, or is treated by a secondary 
analysis. In the case of the Heliogyro, chordwise camber requires special 
consideration as will be shown later. Extensional strain parallel to the blade 
axis is neglected on the assumption that the extensional rigidity is very large 
compared to the bending and twisting rigidities. Vertical and inplane shearing 
flexibilities are similarly ignored. 


All of the above assumptions, except as noted, self-evidently apply to 
the sample Heliogyro blade whose properties are listed in Table 2. In addition 
the thinness of the blade permits us to ignore the elastic bending stiffness in 
a plane normal to the blade chord and the elastic torsional stiffness, in com- 
parison with the corresponding centrifugal stiffening effects. 


An uncoupled analysis of blade motions is one in which the differential 
equations for twist, 0, vertical deflection, w, and inplane deflection, u, are 
uncoupled. In such an analysis the coupling due to off-axis center of gravity 
location, and the nonlinear effects of finite motions, etc., are ignored. The 
uncoupled equations yield valuable information regarding vibration modes and 
the dynamic response to excitation. The uncoupled small-motion equations, 
which incorporate the assumptions which have been mentioned, are: 


TWIST: 


© © © 

- hfrr If) + e ^ 19 * pI If ■ "e 


( 6 ) 
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VERTICAL DISPLACEMENT: 



The symbols that occur in the above equations are defined in the list of 
symbols. In particular, a r is the spanwise component of stress due to centrifugal 
force. The symbol s(7), (F), (3) , and(F) classify the terms according to their 
origin as fol lows: 

© - terms derived from elastic stiffness 

© - spanwise stiffness terms derived from rotation (centrifugal 
st i f fen i ng) 

(?) “ local stiffness terms derived from rotation (gyroscopic stiffness 
terms) 

Q) - ordin ary inertia terms.. 

The terms on the right hand side represent the loading per unit spanwise 
length and also the coupling effects that will be developed later. It will be 
noted that only the inplane displacement contains an elastic stiffness term. 

The gyroscopic stiffness term in the first equation (which is sometimes called 
the tennis racket effect) is stabilizing in the sense that it resists applied 
moment. The gyroscopic stiffness term in the inplane displacement equation is 
destabilizing but the combined effect of the second and third terms is stabilizing 
(or neutral in the case of rigid rotation about the z axis). The equations have 
been written so that they are correct when the blade properties (p, h, c, I, E) 
vary along the span. The blade properties will, however, be assumed to be con- 
stant in most calculations. 

A set of boundary conditions is required in addition to the differential 
equations. At the outboard end the boundary conditions are that no concentrated 
forces are applied. At the inboard end the boundary conditions depend on the 
nature of the body located at the center, and on the number of blades connected 
to it. It is convenient, for initial treatment, to assume that the center body 
is massive so that translations (u, w) and the inplane rotation, £, are zero. 
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If many blades are present, these conditions are also correct in the presence of 
finite hub mass for many of the vibration modes (all except those with so-called 
cyclic and collective symmetry, ref. 10). The root end boundary condition on 
the pitch angle, 6, is that it must follow the motion imposed on i t by the pitch 
control mechanism. 

The boundary conditions should be applied at a short distance from the 
center line due to the finite length of the attachment fittings. Such offsets 
are important for helicopter rotor analysis. Due to the extreme length of 
Heliogyro blades, the boundary condition offsets are only a very small fraction of 
blade span. Their effects will be treated by perturbation techniques applied 
to the solutions for zero offset. 


It is convenient, for the study of trends, to reduce the equations of 
motion to dimensionless form. This may be done by introducing the following 
scaling for the independent and dependent variables. 


Independent variables: 



(9) 


T = Qt 


( 10 ) 


Dependent variables: 



( 11 ) 


(12) 


In addition it will be assumed that the blade properties are uniform along the 
span. For that condition the spanwise stress due to centrifugal force is 

a r = 1 pft 2 (R 2 - r 2 ) (13) 


Dimensionless equations of motion are then obtained by dividing equation (6) by 
pfi 2 I and equations (7) and (8) by pf2 2 chR. The results are: 


TWIST: 



+ 0 + 


3 2 9 

3T2 


O'!) 
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VERTICAL DISPLACEMENT: 


2 9y V Y 9y 


+ lin = 7 

3t2 t z 


05 ) 


INPLANE DISPLACEMENT: 


1 3 " 


2 3y2 V 3y2 


3 2 u 


1 !_((i- y 2)3u. u + ify. = 7 

2 3y I 3y / 9T 2 x 


( 16 ) 


where m_, f , and f are the scaled loads and K is a dimensionless parameter, 
0 Z X 


= El 
j pft 2 chR 4 


_L JL 

12 O 

o 



( 17 ) 


is the spanwise stress at the blade root due to centrifugal force, 


a 


o 


= y pft 2 R 2 


( 18 ) 


K, which is the only parameter appearing on the left hand side of the 
equations of motion, is a measure of the elastic chordwise bending stiffness 
relative to the stiffness provided by centrifugal force. Since equation (16) 
is dimensionless and since the range of y is 0 < y < 1, the magnitude of K 
provides a direct indication of the relative magnitude of the elastic bending 
stiffness. For the sample blade whose properties are listed in Table 2, 

-4w'(wT-- 8 »^o-‘ OS) 

For comparison, the value of K for a conventional helicopter blade is of the 
order of unity. 
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UNCOUPLED VIBRATION MODES 


The homogeneous equations of motion, when written in dimensionless form, 
contain only a single parameter, and it has been shown that that parameter, K, is 
small for Heliogyro blades. If K is set equal to zero and if sinusoidal motion, 
e' WT , is assumed, the uncoupled homogeneous equations become 


TWIST: 

+ e = 0 (20) 

VERTICAL DISPLACEMENT: 

537 ( 21 ) 

INPLANE DISPLACEMENT: 

537( (, -* z C) + ^ * ') 0 < 22 > 

All three equations are variations of Legendre's equation, for which the 
standard form is 


d_ 

dy 


(< 


_y2 ^) + n ( n+, ) x = 0 


(23) 


The parameter n in equation (23) may be considered to be an eigenvalue whose 
value is determined by boundary conditions. For the range, 0 < y < 1, the eigen- 
values for zero displacement at y = 0 are the set of odd integers, and the eigen- 
values for zero slope at y = 0 are the set of even integers. The corresponding 
eigenvectors are the Legendre polynomials, (y) , listed below. 


n 

p „M 

P n <°) 

p n m 

0 

1 

1 

1 

1 

y 

0 

1 

2 

j (3y 2 - 1) 

-1/2 

1 

3 

j (5y 3 - 3y) 

0 

1 

A 

3- (35y 4 - 30y 2 + 3) 

3/8 

1 

5 

g- (63y s - 70y 3 + 15y) 

0 

1 
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Formulas for the vibration mode frequencies, obtained by comparing equation 
(23) to equations (20), (21) and (22), are: 


TWIST: 


U) 


i± = -> / n(n-H) ~ 

Q V 2 


(2A) 


VERTICAL DISPLACEMENT: 


— _ W / n (n+1 ) 

w t i V “ 


(25) 


INPLANE DISPLACEMENT: 


to 


u> 

a 



n (n+ 1 ) , 

2 


( 26 ) 


The scaled frequency, to, is the value of the vibration frequency in cycles 
per revolution. Table 3 gives the numerical values of oj for the first four modes 
of each type with zero root end motion. 

It is noteworthy that the vibration mode frequencies bear a fixed relation- 
ship to rotor speed, regardless of size, shape or blade material properties. It 
is fortunate, from the viewpoint of resonant excitation due to environmental 
effects, that none of the higher (n > 1) vertical and inplane mode frequencies 
are exact harmonics of rotor speed. The twist mode resonance at ^/rev. is not 
serious because damping can easily be introduced into the pitch degree of freedom, 
as will be shown . 

The lowest vertical and inplane modes at n = 1 correspond to rigid body 
rotations of the blade about a fixed hub. The vertical mode frequency, which is 
at 1/rev. in a coordinate system that rotates with the blade, occurs at zero 
frequency and also at 2/rev. in a nonrotating coordinate system. These frequencies 
correspond to the precession and nutation of the vehicle as a whole. 
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Only two of the solutions with even order Legendre polynomials have 
practical significance. They are the lowest order (n = 0) twist mode which has 
a frequency of one cycle per revolution, and the lowest order vertical mode 
whose frequency is zero. The vertical mode is trivial since it simply describes 
rigid body translation parallel to the axis of rotation. The existence of the * 
twist mode, on the other hand, demonstrates the important result that a cyclic 
pitch imposed at the blade root at a frequency of one cycle per revolution will • 
result in a uniform response along the blade span. This kind of excitation is 
employed in several types of maneuvers, see reference 9- 

At excitation frequencies other than 1/rev., the pitch response to excitation 
imposed at the root is not uniform due to the imbalance between the inertia and 
the local centrifugal stiffening terms. Figure 5 is a plot of the ratio of blade 
tip response to blade root excitation as a function of the exciting frequency, 
obtained by numerical integration of equation (20). The result for zero excita- 
tion frequency (0 . /0 = .319) is important because it is the steady state 

^ 'tip root 

limit for collective pitch excitation. 

The solutions that have been presented for uncoupled inplane modes are sub- 
ject to the assumption that the elastic chordwise bending stiffness is zero. 

Except for the lowest mode, the error due to this assumption is unimportant for 
the very slender blades envisaged for the Heliogyro. The error in the lowest 
(n = 1) mode depends on the boundary condition for chordwise slope at the blade 
root. For modes with collective symmetry (modes in which all blades execute 
identical motions) the error is zero because the appropriate boundary condition, 
which is that the chordwise bending moment is zero, is satisfied by the assumed 
eigenvector. For all other types of modes, the appropriate boundary condition is 
that the chordwise slope is zero, which results in a finite value for the lowest 
mode frequency. 

An analytical expression will be derived for the lowest mode frequency on 
the assumption that the chordwise bending stiffness parameter, K, is small but 
finite. The approach will be to derive an expression for the chordwise deflection 
that satisfies the boundary conditions at the root (u = 0, du/dy = 0) , that 
approximately satisfies the differential equation near the root, and that 
asymtot i cal ly approaches the straight line (u = y) predicted for K = 0. The 
expression for deflection will then be substituted into a Rayleigh quotient to 
obtain an approximate frequency. 
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Neglecting inertia and local centrifugal stiffening terms in equation (16), 
the homogeneous differential equation for inplane bending near the root, y « 1 , 
i-s approximately 


K 


d_jj 
dy ^ 



(27) 


A solution which satisfies the differential equation and the boundary 
conditions, and that remains finite as y increases is 



Note that, for y » , the deflection is approximately a straight line 

that intersects the y axis at y = yv. 

The kinetic and potential energy functions for the blade are easily identi- 
fied by examination of the terms in equation (16). They are 



The Rayleigh's quotient is 



o 
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If equation (28) is substituted into equation (30, 
as a small parameter such that higher order terms may be 
with lower order terms, the resulting expression for the 
i s 


and i f ST is treated 
neglected in comparison 
lowest inplane frequency. 





(32) 


Straightforward analysis, reference 11, of a uniform articulated rotor blade 
with lag hinge offset, e, shows that the lowest inplane frequency is (approximately) 



By comparison of equations (32) and (33), the "effective" lag hinge offset 
for a Heliogyro blade is 


e = R/iT ( 3*0 

For the Sample Blade design, K= .833 x 10 **, so that co = ^( . 833x1 0 4 ) ^ 

= .1182/rev. From the form of equations (17) and (32), it is seen that the 
frequency of the lowest inplane mode is inversely proportional to the square root 
of the blade aspect ratio, R/c. 

A rough estimate of the effect of chordwise bending stiffness on the higher 
modes can be obtained by substituting the assumed displacement function for K = 0 
into the potential energy function, equation (30). The resulting frequency is 


co = to Q (l + e) 


(35) 


where co q is the frequency for K = 0 and 


2v c - 2 r-T 


-H00 


■f G? 


- 2u 2 


( 36 ) 


dy 


As an example consider the second mode (n = 3) for which the deflection function 
(Legendre polynomial) is u = 5y 3 - 3y. 
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For this function, 


and the result is 


2 


(u) = 25y G - 30y“ + 9y 2 


(^) = 225y‘* ~ 90y 2 + 9 


(— ) 

\dy 2 / 


900 y‘ 


( 37 ) 


£ = 52. 5K 


(38) 


For the Sample Blade, K - .833 x 10 and e = .00436. The importance of the 
chordwise bending stiffness increases as the mode index increases. For high 
order modes the eigenfunction is approximately a sine wave, 


u = s in 


2Tty 


(39) 


where X is the wavelength. An estimate of the effect on frequency is obtained 
by neglecting 2u 2 and replacing 1-y 2 by 2/3 in the denominator of equation (38). 
Thus 


e 


'(F 



' 2ir\‘ 

3 1 

a ) 


3tt 2 K 

X 2 


(40) 


The condition under which the elastic energy and the potential energy of centri 
fugal stiffening are equal, is £ = V e /2V = 1/2, or 


X = -n/SK 7 


(41) 


— 4 

For the Sample Blade, K ** .833 x 10 and X = .0722. The number of nodal points 
along the span is approximately equal to 2/X. This is a large number for the 
Sample Blade, again indicating that the elastic bending stiffness has a negligible 
effect on frequency in the practical range of excitation. 
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DAMPING DEVICES 


Heliogyro blades have almost no inherent damping. Aerodynamic forces, 
which provide most of the damping for conventional helicopter rotors, are 
totally absent. The effectiveness of internal structural damping is small 
because elastic stiffness accounts for only a small part of the total stiffness. 
In fact, only the inplane modes have appreciable structural damping and then 
only in the higher modes. Damping for all other blade modes must be supplied 
by damping devices. 

Damping is required in structures primarily for two reasons. The first is 
to limit the amplitude of the response of the vibration modes to sustained tran- 
sient excitation (which over a long period of time may be regarded as either a 
periodic or as a random excitation). The second is to suppress instabilities 
involving structural deformation which may result from control feedback or non- 
conservative environmental processes. 

Both of the above reasons apply to the Heliogyro. A primary source of 
transient excitation is the blade pitch control mechanism which can execute 
abrupt motions in response to maneuver commands. Another source of dynamic 
excitation is the solar illumination which can, in conjunction with cyclic 
blade pitch, produce periodic radiation pressure fluctuations at several of 
the lowest harmonic orders of rotor speed. Fortunately, as has been shown, 
there are no prominent resonances of vibration modes with the lower orders of 
rotor speed. Gravity gradient also produces periodic load fluctuations of 
relatively small size. 

The primary means for combatting periodic excitation is to avoid resonances 
between exciting frequencies and natural modes. As long as this is done, 
periodic excitation produces no requirement for damping. The response to random 
excitation, on the other hand, is inversely proportional to the square root of 
the damping (as measured by the ratio of the modal damping force to the modal 
inertia force). Thus, even a very small amount of damping is effective in re- 
ducing the amplitude of response to sustained transient excitation. 
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The Heliogyro is potentially subject to structure-autopilot instability 
because it requires an autopilot in order to maintain trimmed flight and to 
e.xecute coordinated maneuvers. Even though the gain of the autopilot may be 
very small at frequencies corresponding to the higher vibration modes, some 
small level of structural damping is required in order to avoid instability. 


Thermal expansion is a source of instability that has been troublesome for 
long booms on space vehicles. This mode of instability involves differential 
thermal expansion between the illuminated and non i 1 1 umi nated sides of the boom, 
which produces a change in the illumination incidence angle which, in turn, 
changes the temperature distribution. Thermal time lag plays an important role 
in determining the phase of the feedback. Fortunately for the Heliogyro, 
differential thermal expansion between the leading and trailing edges of the 
blades is very small so long as the blades are reasonably flat, and the thermal 
time lag is small due to the thinness of the blades. Nevertheless, some small 
amount of damping is probably required to prevent this kind of instability. 


Since it has been determined that damping devices are required, it is 
pertinent to seek practical ways of providing them and to estimate the damping 
levels that can be achieved. The most desirable locations for the damping 
devices are at the ends of the blade. Distribution of devices along the blade 
is undesirable because they will increase the bulk of the stowed configuration 
and will reduce the reliability of deployment. The blade tip is not a practical 
location for damping devices because significant forces cannot be produced there 
except at the expense of a significant tip mass. The most practical location 
for damping devices is at the blade root attachment, which, incidentally, is the 
location chosen for conventional helicopter rotor blades. 


As long as the motions of the blade are substantially uncoupled, separate 
devices must be provided for pitch, vertical displacement and inplane displace- 
ment. A higher level of damping can be provided for pitch than for the other 
degrees of freedom, which is fortunate because of the important role played by 
pitch motions in maneuver control. 
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Pitch Damping 


Damping is conveniently provided for the pitch degree of freedom by means 
of a parallel combination of a spring and a dashpot between the control actuator 
and the blade as shown in Figure 6. In order to achieve maximum damping, the 
damper spring should be of the order of magnitude of the blade torsional restraint 
due to centrifugal force (pf2 z lR), and the time constant of the damper (B^/Kj) 
should be of the order of the period of one of the lower torsional modes. A 
parametric investigation was made of the idealized model shown in Figure 6 in 
order to determine an optimum damper combination. The criterion used was the 
number of revolutions for the blade tip motion to damp to within 10% of its 
steady state value in response to a collective pitch input. The collective 
pitch input was a ramp function that achieved its final value in one-half 
revolution, as shown in Figure 6. Results for a typical case are shown in 
Figure 7 where the blade tip response, 0^^, 5r, d a l sc the blade root response 
just outboard of the damper, 0^, are plotted versus time measured in rotor 
revolutions. The results for all cases are summarized in Figure 8 where it is 
seen that the selection of the damper parameters is not critical. 

An eigenvalue analysis was also made made of the configuration whose 
transient results are shown in Figure 7 in order to determine the effects of 
the damper on the torsional vibration modes. The results are presented in 
Table A. 


The table shows, as does Figure 7, that only the lowest mode has significant 
response after two revolutions. Comparison of the frequencies listed in Table A 
with the results for the blade without damper given in Table 3> shows that the 
damper slightly reduces the frequencies of the lowest two modes. The increases 
shown for the 3rd and Ath modes are the result of finite element approximation. 

It is instructive to calculate the physical magnitudes for the damper and 
the spring for the Sample Blade whose properties are listed in Table 2. The 
control stiffness is 

K , = 2pft 2 IR = z- ft 2 pc 3 hR 
d 0 


1 

= T CT 

3 o 


cfh 

R 


= JX (1000) 


(120) 3 x .25 x 10" 3 
120,000 


m 


1.2 in-lb/rad. 
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The control damping is 


B d ■ t K d ’ :1m x K2 ’ 242 ln - ,b/rad/sec < 43) 

The control stiffness, and consequently also the control moment, is extra- 
ordinarily small in view of the size and weight of the rotor blade, indicating 
that the pitch control mechanism will be delicate and that friction about the 
control axis must be minimized. Perhaps the best way to achieve the required 
control stiffness and damping is by means of the output impedance of an electric 
motor which can, with appropriate servo control, be adjusted to virtually any 
value. The effect of friction can be minimized simply by increasing the magnitude 
of the control stiffness at the expense of a decrease in the damping decay rate, 
as shown by Figure 8. A design concept for a complete Heliogyro blade retention 
system will be introduced after the means for damping inplane and vertical blade 
motions have been discussed. 


Inplane Damping 

Inplane bending deflections are the only blade motions that have appreciable 
structural damping. The theory presented earlier for the effect of elastic 
stiffness on inplane mode frequencies can also be used to estimate mode damping. 
If g g is the structural damping coefficient of the elastic material, then the 
modal damping parameter is, approximately, 

9 = /g e (M) 

c 

where V /V is the ratio of the elastic energy to the potential energy of centri- 
fugal stiffening. From equations (36) and (38), the energy ratio for the second 
inplane mode is 



(45) 


where K is the inplane bending stiffness parameter. For the Sample Blade, 
K = .833 x 10 \ so that 
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(W) 


g = 105 x .833 x 10"^ g e 

= .00872 g g 

Thus, the effective damping of the second mode is less than 1% of the 
material damping. The effectiveness of structural damping increases rapidly 
for higher modes, as indicated by equation (40). It is also high for the lowest 
mode but it cannot be considered adequate due to the very low natural frequency 
of this mode. 

A practical design concept for an inplane damping device is a damper that 
resists relative rotational motion between the blade root and the center body. 
The damper must have some stiffness in order to maintain the azimuth angles 
between blades (provided that there are more than three blades). 

An idealized model for the inplane motion of the rotor blade when vibrating 
in its first mode is shown below. 




Damper 

■"TRRRT'- 

K d (l+is d ) 


o 

-o- 


Blade 


-nm ^ — o 


The damper is represented by a spring whose value is given an imaginary 
part to represent hysteresis damping. The blade is represented by its inertia 
for rigid body rotation about the blade root, 1^, and by a spring K^, that is 
calculated from the inertia and the frequency of the lowest mode 


K c - U *I C 


( 47 ) 


Using dimensionless coefficients, 


1/3 


( 48 ) 
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and, from equation ( 32 ) 


a) 2 = I /T 7 (49) 

Thus, in dimensionless form, 

\ = J ^ (50) 


The damping of the first mode is, approximately, for g^ « 1, 

9 i V , + V 9 d ( 5 

a C 

where is the energy stored in and is energy stored in K^. Since the 
energy distribution of springs in series is inversely proportional to their 
st i ffnesses 


9 i K . + K 


(52) 


Evidently maximum damping is achieved when is small compared to K^. 

An idealized model for vibration in the second inplane mode is shown below 



The blade is represented by its modal mass, I^, and modal stiffness, K 2> 

referred to the slope, r , at a point near the root. The equivalent elastic 

b 

coupling spring, K , that was calculated from the properties of the first mode, 
is connected between the damper and the blade. The modal mass is calculated 
from the deflection function u = £^(y - j y 3 ) and the kinetic energy equivalence, 


T = 


2 





(53) 
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Thus, in dimensionless units, 


= / (y "f y3)dy =fe 


(5*0 


and 


„ _ -2, c w ** _ 20 

K 2 oj 2 I 2 5 x 


(55) 


The damping of the second mode is, for g^ « 1, 

V d 

9 2 V , + V + V 0 9 d 
d S 2 


(56) 


where, in addition to previously defined terms, V 2 is the energy stored in K 9 . 


Let 


K = 


K ^ K d 


S K c + K d 


2 ' 

(57) 


be 


the equivalent stiffness of K, and K in series. Then 

d £ 


V d + \ + U 2 - T ^< K 2 * K C > 


and 


(58) 


V + v = — r 2 K 
d \ 2 ‘’b 


so that, using equations (51), (52) and (58) 


= 


V , + V 
d _L_ . 


2 V . + V V . + V + V. 
d £ d £ 2 

K 1< 

= K . + K * K + X * 9 d 
d C 2 C 


(59) 


_L 


Kd + K? , + K 4 + f 
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The maximum value of obtained by differentiating equation (59) with respect 
to K^, occurs when 


K 


d 



( 60 ) 


For very slender blades, « 
damping is obtained by setting = 
and (59), 


l^, so that nearly optimum second mode 
K . Using this value in equations (52) 


9 1 2 9 d 


9 2 = 


2 9 d 


2 9 d 


2K_ 


1 + 


1 + 


£30 


For the Sample Blade /K" 7 = .00931 , so that 


9 2 = 


2 9 d 
137 


( 61 ) 


( 62 ) 


A typical value of g^ for a good elastomeric material, such as silicone 
rubber, is 0 . 30 . For this value of material damping the modal dampings for the 
Sample Blade are g^ = .15 and g^ = . 00109 . The corresponding times to half 
amplitude, measured in revolutions, are 


t 


1 


1 n2 
iTg ] aj 1 


12.5 revs 


t 


2 


1 n2 

TTg 2 ^2 


= 90.7 revs 


(63) 


The relatively long decay time for the 2nd mode is a direct result of 
blade slenderness as expressed by the elastic stiffness parameter /iT" 7 . From 
the definition of K, equation (17), 


/r-fV 



o 


(61.) 
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Thus, other factors remaining equal, the damping of the second inplane mode (and 
of higher modes) is approximately proportional to the ratio of blade chord to 
blade span. The inability to provide adequate damping for inplane modes may set 
a limit to blade aspect ratio. 

Somewhat better results can be achieved with viscous dampers. A configura- 
tion using two springs and two viscous dampers is shown below 


dl 


d2 







K dl 

\AA/ 


K d2 

V\A/ — 1 


An eigenvalue analysis was made for the inplane vibration modes of the 
Sample Blade design using this damper configuration with the following dimen- 
sionless component values 


K 


dl 


K d2 = 2K C 


“l B dl 


dl 


k 

3 


U) 2 B d2 


d2 


1 

2 


(65) 


The springs were selected to be equal and to give the same stiffness as 
the hysteresis damper described earlier. The dashpots were selected to provide 
maximum damping for the first and second modes. The results of the analysis are 
shown in Table 5. It will be noted that the presence of the damper slightly 
reduces the lowest mode frequency (see discussion following equation (3^)). The 
increase in the higher mode frequencies over the values for the blade without 
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bending stiffness, see Table 3, is partly due to bending stiffness but mostly 
due to finite difference approximation. The damping for the first mode is 
slightly less than that for the silicone rubber damper but the damping of the 
second mode is substantially larger. Unfortunately the vibration frequencies 
for real blades are so low (their periods are of the order of minutes) that 
the physical realization of viscous dampers may not be practical. 

As in the case of the pitching modes, it is instructive to estimate the 
physical magnitudes of parameter values. The physical value of the damper 
spring for the Sample Blade is 


K* = pchft 2 R 3 K d 

-K)W(? >*’) 

= 2000 X 36 OO x.00^655 

= 33600 in-lb/rad. ( 66 ) 


This stiffness can, for example, be provided by an extensional spring with 
a spring constant equal to 336 lb/in located ten inches from the hinge axis. 

I t wi 1 1 be shown later that a damper containing less than one cubic inch of 
silicone rubber will satisfy strength requirements. 

Vertical Damping 

A problem encountered in designing a damper for vertical blade motions is 
that the blade has essentially no bending stiffness to react the moment applied 
by a rotational damper. Effective bending stiffness can, however, be provided 
by structural reinforcement over a short span of blade at the root end. A 
design concept is shown in the following sketch. 
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DAMPER 


3 RACE 



The reinforcement consists of a brace made from two thin sheets of material 
similar to the blade material. Additional spanwise reinforcing fibers are 
attached to the sheets. In the stowed configuration, the bottom sheet is back- 
rolled under the top sheet. 

If the brace and the damper were infinitely rigid, the blade would be 
provided with a flap hinge offset at a distance r b from the roller, and the 
flapping moment at the roller would be 


M 

r 


r, T 3 
boo 


(67) 


where T q is the centrifugal tension and is the flapping angle just outboard 
of the flap hinge. An idealized model for vertical blade motion, including an 
elastomeric damper and a brace with finite stiffness is shown below. 


DAMPER 


-Wv 

K d (l+ig) 







( 68 ) 
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In the model I. and K. are the modal inertia and stiffness of the i th mode 
I I 

W h 

referred to 3^ = — , the normalized displacement at the hinge. K is the 

b 

centrifugal stiffness associated with the increment of span occupied by the 
brace. Its value, in physical units, is 


K" = r, T 
c bo 


(69) 


or, in dimensionless form, 


K c 2 


(70) 


where y^ = r^/R. The value of the brace stiffness is, in physical units 


(EA) b d 2 

2r7 


(71) 


and, in dimensionless form 


(EA) 


K b = 




a ch \ R / 


(72) 


(EA)^ is the extensional stiffness of one sheet. 

The modal inertia, I., can, as in the case of inplane bending, be evaluated 
from the spanwise integral of kinetic energy. Assuming that the mode shapes are 
adequately described by the Legendre polynomials discussed earlier, the dimen- 
sionless values of the modal inertias and stiffnesses for the first two modes, 
referred to blade flapping, 3 q . just outboard of the hinge are 


10 


X 20 


K 10 = W 1 2I 10 


( 73 ) 
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W 2 2l 20 


2b 
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A transformation of these quantities to the normalized hinge displacement, 
is provided by the ratio of the moment on the roller, M^, given by equation 
(67) to the generalized force acting on 3 o> 


M 

o 


= K. 3 
1 00 


(7*0 


In dimensionless form, the moment on the roller is 



(75) 


The blade Inertia, I., and stiffness, K. , in the idealized model are, there- 
fore, 



(76) 


Substituting from equation (73) into equation ( 76 ) 


1 1 = i y b 2 


I = air y 2 
2 57^ Y b 


K 1 " i y b 2 


(77) 


K 


2 



2 


As in the case of inplane motion, the ratio of modal damping to material 
damping is approximately equal to the ratio of the energy stored in the damper 
spring to the total potential energy. Thus 


Vd 


V. + V. + V + V. 
d b c 1 


(78) 
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The simple idealized model can readily be analyzed to produce the energy 
ratio in terms of stiffnesses. If it is assumed that y^ « 1, then, to a good 
approximation, 


9i ' 9 < 


(79) 


where = ^K^/K^K^. Cursory examination of equation (79) shows that maximum 
damping is obtained by making large. This may result in a significant weight 
penalty so that selection of an optimum design involves a tradeoff between weight 
and damping. 


As an example for the purpose of estimating the order of magnitude of 
damping that can be expected, consider the Sample Blade with a brace described 
by the following parameters: 


d = 1.0 ft 


r b = 50 ft 

= 30 x 10 6 psi (steel) 
A, = .010 in 2 

D 


Then the dimensionless stiffness of the brace is 


K b = 


itx.005 


300,000 

30 


-faY ■ 


50 x 10 


-k 


The weight of the brace is 3-6 lbs which is 1.9% of the blade weight. The 
blade stiffness parameters are, from equations ( 70 ) and (77), 

Kj = . 1875 x lO' 4 

K 2 = ,161»2 x 10"' 4 
K. = 25 x 10 _it 

An optimum value of can be found for given values of and K c< For the 
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parameters that have been assumed, the optimum value is near 


K , = 20 x 10 
d 


The corresponding values of modal damping are, assuming an elastomeric 
damper with g^ = 0 . 30 , 

g, = 3.73 x 10 _Zt 
g 2 = 3.27 x lO' 2 * 


These values are only 
inplane mode. They can be 
length and brace weight or 
diameter of the roller can 
extension of arms that are 


about one-third as large as the damping for the second 
increased at the expense of an increase in brace 
by an increase in the diameter of the roller. The 
be increased, in the deployed configuration, by the 
stowed inside the roller. 


The physical value of the damper stiffness is 


Kj = 1 A , 2»00 in-lb/rad 


which is about four-tenths 


of the value of the inplane damper stiffness. 


Damping of Modes that Involve 
Rotations of the Center Body 


If all blades are identical, certain types of modes will exist in which 
the inplane and vertical dampers discussed above are ineffective. They are 
modes in which, due to the symmetry of the blade motions, the blade bending 
moments are zero at the center of rotation. There are two classes of such modes 
as shown in Figure 9. 

1 1 wi 1 1 be noted that the modes of each class involve a rotation of the 
center body. It would, therefore, appear feasible to damp them by attaching 
damping devices to the center body, or by attaching devices, such as gyroscopes, 
that provide a significant impedance against which the blade dampers can react. 
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Unfortunately the size of the center body will probably be so small compared 
to the blade radius, that appreciable damping cannot be provided in this manner 
without a large weight penalty. 

A more fruitful approach is to alter the characteristics of the blades so 
that the bending moments at the hub do not vanish in any modes, except, of 
course, in the lowest modes which correspond to rigid body motion of the entire 
vehicle. For a two-bladed rotor, this can be accomplished by making one blade 
longer than the other so that the center body is no longer at the center of 
rotation. For multibladed rotors, a convenient technique is to taper the chord 
of some of the blades so as to alter their mode frequencies- (Simply changing 
blade length will not significantly alter the mode frequencies-) For example, 
the second vertical mode of a blade whose chord is linearly tapered to a value 
of zero at the tip is 1.955/rev compared to 2.^49/rev for a uniform blade, and 
the second inplane mode is 1.68/rev compared to 2.236/rev. These differences 
are large enough so that, considering the relatively small influence that hub 
rotational constraints have on blade mode frequencies, the modes of the tapered 
and uniform blades will be virtually uncoupled. Thus, one set of blades will 
be nearly motionless in the modes of the other set and the dampers will be 
effective. An attractive configuration is a symmetrical six bladed rotor with 
three uniform blades and three tapered blades. Rotors with fewer blades must 
be unsymmetri ca 1 in order to make the dampers effecti ve for all modes. 

One last mode of interest which is not damped by any of the techniques 
discussed so far, is the nutation mode of a vehicle with three or more blades, 
in which the entire vehicle vibrates as a rigid body at approximately two 
cycles per revolution. In the rotating system this mode appears as a blade 
flapping mode at one cycle per revolution. Active control feedback can be 
used to damp the nutation mode. Also it may be feasible to damp the nutation 
mode by means of a tuned absorber that consists of a long boom attached to the 
center body in the direction of the spin axis. 

A solution to the damping problem which has the advantage that it applies 
to all modes that involve rotations of the center body, including the nutation 
mode, is to employ an auxilliary set of counter-rotating blades. The counter- 
rotating blades, which must have their own dampers, need only be large enough to 
provide an appreciable impedance to hub rotation. They can be fabricated from 
high strength materials in order to minimize weight. 
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BLADE RETENTION SYSTEM 


The need for blade dampers about all three rotational axes has been 
identified and appropriate values of stiffness and damping have been calculated 
for the Sample Blade design in the preceding section. The recommended properties 
of the dampers are as follows: 


Mot i on 

Damper Type 

Spring 
Constant 
( in- 1 b/rad) 

Damping Constant 

P i tch 

Electric Motor 

1.2 

B, = 2*t2 i n- 1 b/rad/sec 

d 

Vert i ca 1 

S i 1 i cone Rubber 

1M00 

cn 

CL 

II 

0 

Va) 

1 np 1 ane 

S i 1 i cone Rubber 

33,600 

g d = 0,3 


It is evident, from comparison of the spring constants that the friction 
about the pitch axis should be as small as possible. Additional requirements 
for the blade retention system are that it support the centrifugal load (30 lbs), 
that it accommodate blade motions of the required amounts, and that it withstand 
environmental hazards. The major environmental hazards are shock and vibration 
(during launch) and micrometeoroids. 

The available pitch motion should be large in order to maximize maneuver 
capability. At least A5° should be provided. The required vertical flapping 
accommodation can be estimated to be a reasonable factor times the steady coning. 
Table 2 records the blade root tension ( 30 lbs) and the total thrust for normal 
solar illumination (.0188 lbs) on the Sample Blade. Thus, assuming that the 
center body is much heavier than the blades, the root coning angle is 

8 = = .000626 rad 

o 30 

It seems reasonably conservative to assume dynamic flapping of equal magnitude 
so that the total required flap angle accommodation is approximately .0013 radian 

= .075°. 
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The required inplane bending moment can be estimated from the analysis of 
maneuvers. A reasonably severe maneuver (sudden application of cyclic pitch) is 

studied in a later section of this report. The results show that the peak in- 
plane bending moment is 13.5 in-lbs per radian of applied cyclic pitch. It seems 
reasonably conservative, therefore, to take twice this value, 27.0 in-lbs, as 
the design inplane bending moment for the damper. The corresponding angular 
deflection of the damper is 


AC 


27 

d 33,600 


.000802 


rad . 


We are now able to calculate 
the vertical and inplane dampers. 

V = ■ 


the required volume of silicone rubber for 
From the equivalence of two energy expressions 


ll 2 Vol = lUi- 
2 G 2 K 

d 


(80) 


the required volume of material is 


Vol 



( 81 ) 


where G is the shear modulus and x is the allowable shear stress. The shear 
modulus for silicone rubber is approximately 100 psi. A conservative estimate 
of the fatigue endurance limit is 5 psi. The required volume of material for 
the inplane damper is, therefore, 


Vol 


100 (27. 0) 2 

I5P" X 33,600 


.087 in 3 . 


The flapping moment at the roller for a rigid damper and a rigid brace is, 
from equation (67), 


M r = r b T Q 8 o = 600 x 30 x.0013 = 23. in-lb. 

For the damper design studied in the previous section, the moment applied to the 
flexible damper is 36% of this value, or 8.** in-lb. The required volume of 
material for the vertical damper is, therefore. 
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Vo ' ■ - 020 in ‘- 

A design concept for the blade retention system is shown in Figure 10. 

It is similar, in many respects, to retention systems currently in use on 
several light helicopters. Centrifugal tension is taken by three stranded wire 
cables which can accomodate large pitch motions without appreciable torque. The 
maximum pitch angle is approximately ±55°* The assembly includes only a single 
bearing which carries no thrust. The values of the design shear loads for the 
bearing, which are equal to the design bending moments for the damper divided by 
the moment arm (ten inches), are 2.7 lbs. inplane and .84 lbs. vertically. 
Assuming that the bearing radius is 1/4" and that the coefficient of friction is 
0.02, the maximum friction torque is 

T = *Y(2.7) 2 + (.8^) zl x i x .02 = .014 in-lb. 


The corresponding break-away pitch angle is .014/1.2 = .0116 rad = 0.67°, which 
is certainly tolerable. 

The pitch drive shaft passes through a central hole in the pitch bearing 
and is attached directly to the rotor of an electrical motor. The friction 
loads on the drive shaft should be negligible. 

The pitch case, the meteoroid shield and the droop stop are designed to pro- 
tect the mechanism from environmental hazards. 
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FORCES DUE TO SOLAR ILLUMINATION 


Figure 11 shows the coordinate geometry for a rotating blade that is coned 


through an angle 0 with respect to a tangent to the cone of rotation. A further 
rotation, £, about the normal to the blade produces no change in the radiation 
pressure and need not be considered. The illumination lies in the x,z plane 
at an angle y with respect to the axis of rotation. 

The following expression is derived in ref. 9 for the normal pressure on 
the blade surface under the assumption that the coefficient of reflectivity is 
unity. 

P n = p Q [(sin 0 sin \p - cos 0 sin 3 cos iJj) sin y + cos 0 cos g cos y] 2 ( 82 ) 

where \p = fit, is the azimuth angle relative to the. plane of illumination. The 
components of vertical and inplane load for use in structural dynamics, equations 
(7) and (8), are 


The moment about the pitch axis due to solar illumination is zero provided 
that the reflectivity of the surface is uniform and that the blade is not cambered 
in the chordwise direction. The effects of such nonuniformities will be dis- 
cussed later. 

Equations (82) and ( 83 ) are sufficiently general for most purposes. They 
are, however, nonlinear in the motions (0, g), so that it is desirable for con- 
venience in analysis to derive linearized expressions for the loads under the 
assumptions of small motion. Some expressions are presented in Table 6 for a 
number of different cases. Since g is always very small, case 5 of Table 6 is 
sufficiently general for all practical cases. It shows that, if 0 varies cyclicly 
through a large angle, 0 = 0^ sin \p, and if the illumination is not parallel to 
the axis of rotation, the vertical exciting force will include harmonic components 
with prominent amplitudes up to the fourth harmonic, and that the inplane exciting 
force will include prominent harmonics up to the fifth. 


through an angle g with respect to the plane of rotation and subsequently pitched 



(83) 


f = p c sin' 0 
rx K n 


Another important fact is that, since 3 is very small and may be neglected 
in comparison with 9, the significant coupling is from the pitch degree of free- 
dom to the vertical and inplane degrees of freedom and not otherwise. Even so, 
the magnitudes of the coupling coefficients are very small when expressed in 
dimensionless form. For example, considering case k in Table 6 


rx 


rx pft 2 chR 


P R 

= 6 cos 2 y • 


6 


2a h = 9 o cos2 Y * f- 


(84) 


where 3 Q is the blade root coning angle under the assumptions that the center 

body is much heavier than the blades and that the illumination is parallel to 

the axis of revolution. For the Sample Blade design the value of 3 is 6.26 x 
-It ° 

10 rad. 
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MECHANICAL COUPLING DUE TO BLADE DEFLECTIONS 


The equations of motion for a rotating beam that is slightly curved include 
small terms that couple the vertical, inplane, and pitching motions. Since the 
curvature of Heliogyro blades is small, the coupling due to curvature would be 
only a slight concern if the torsional stiffness were not also very small. It 
is seen, however, from equations (6) and (8) that the ratio of the effective 
torsional stiffness (o r l) to the inplane bending stiffness (El) is equal to the 
spanwise strain which, of course, is small. Thus, a detailed analysis is required 
in order to determine the significance of the coupling terms. 

It is shown in the Appendix that the pitching moment per unit length due 
to curvature may be expressed as 


m -Mm - 15. M 
0 3r ? 3 r 8 


(85) 


where 8 and £ are vertical and inplane slopes of the blade (see Figure A) and 

M and M are the vertical and inplane bending moments. The bending moments 

P c 

are related to blade curvature via the bending stiffness. Ignoring the flapwise 
bending stiffness about an axis lying in the blade chord, and assuming that the 
pitch angle, 0, is small, the relationships are 


\ - 6 If) + " £ 


( 86 ) 

(87) 


Substitution of equations (86) and (87) into equation (85) produces the 
result 


= Eie fey fey] + Ei|i |fi. + M 

[\3 r / \3r/J 3r 3r £oV3r 3r/ 


( 88 ) 


Equations (86) and (87) are correct to second order in the motions (^, 8> 0) 
and equation (88) is substantially correct to third order. The quantity is 
the chordwise bending moment due to initial stress. Its major source is the lack 
of straightness introduced in the manufacturing process. When the blade is 
deployed it is straightened by centrifugal force, and the initial curvature is 
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converted, almost completely, into a chordwise bending moment. Subsequent pitch- 
ing of the blade rotates the moment to produce a vertical component as well. 

Coriolis forces are another source of coupling between the inplane and 
vertical motions. For an undeflected blade the Coriolis forces acting in the 
inplane and radial directions are 


f x = 2pchfi Tt - 

fy = -2pchfi It 


(89) 


where u and v are respectively the inplane and radial displacements. When the 

9w 

blade is flapped through an angle 8, the velocity normal to the blade, 7—-, has 

~ O L 

dw • . 

a component -3-^- in the radial direction and the radial force, f^, has a com- 
ponent “Bfy normal to the blade. Thus the normal and inplane components of 
Cor iol i s force are 


f cx = _2pch ^ 31 

f c2 - 2„cm |f 


(90) 


The additional normal and inplane loads due to equations (86), ( 87 ) and (90) 


are 


f = 2pchfiB f-H- + - 3 - 
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f x = - 2 P cwg + |L- (Ele - M c 


(91) 

(92) 


The dimensionless forms of the coupling terms given by equations (88), (91) and 
(92) are Q © © 
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(95) 


The significance of each term is determined by comparing its size with the 
terms on the left hand sides of equations (14), (15) and (16). The latter terms 
are of the order of the scaled deformations (0, w, and 7T) . 

The properties of the Sample Blade design will be used in order to obtain 
a quantitative evaluation of the coupling effects. In so doing it will be 
assumed that the vertical deflection is of the same order as the static deflec- 
tion 


3 


s 


*o 

1 + y 


(96) 


and that the inplane deflection is of similar size. The latter assumption will 

be justified later. It will also be assumed that 0 is of the order of unity. 

-h 

The value of for the Sample Blade is 6.26 x 10 rad. Using this value and 
the properties listed in Table 2, it is found that the order of magnitude of 
the terms in equations (93), (94) and (95) are 


term © 
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All of the coupling terms, except possibly those involving M , are small 
compared to the terms in the uncoupled equations, and they may safely be 
neglected. They are, in fact, much smaller than the corresponding terms for a 
conventional helicopter blade. Terms (!) and © are potentially significant 
terms when g and £ have larger values as will be shown later in connection with 
an experimental blade that was tested in a spin chamber. 

As has been mentioned, a major source for M£ q is unavoidable chordwise 
curvature built into the blade during manufacture. Thus, assuming that the 
blade is straightened by centrifugal force during deployment 


1 - « 
V R 


(97) 


where R is the initial radius of curvature. Substitution into the third term 
c 


of equation (93) gives 


s . ■(*»*•») 


(98) 


Static vertical deflection, as given by equation (96), therefore, produces a 
moment load distribution 


03 


• -feXy-A 


(99) 


Figure 12 shows the pitch distribution that results from application of a 

load.m = 1/(1 + y) 2 . Using Figure 12 and equation (99), the pitch angle at the 
0 

tip due to a constant initial radius of curvature is 


0 t = .225 


fexy-. 


( 100 ) 


For the Sample Blade design 


0 = .225 x (500) x 6.26 x 10 


= . 0705 I n~ Jrad ian . 


W 


fe> 
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It should be possible, using only moderately careful calibration procedures, 
to ensure that the initial radius of curvature is greater than 50,000 ft, a 
value which corresponds to a tangent offset of about one inch in 100 ft. Using 
this value 


0 t 


.0705 x 



,0H»1 


rad 


.81° 


A static twist of this magnitude is probably insignificant. Initial curvature 
may be a problem for longer blades, particularly because the coning angle in- 
creases with radius according to the formula 


P K 
60 °o h 


The corresponding pitch angle at the tip is 


0 


t 




( 101 ) 


( 102 ) 


Note that the pitch angle is directly proportional to the square of the radius 
and inversely proportional to the thickness. Note also that it is inversely 
proportional to the elastic energy density, o q 2 /2E, so that materials with high 
energy storage capacity are good candidates for Heliogyro blades. 

Another potential source of steady chordwise bending moment is differential 
thermal expansion. If a linear temperature gradient across the blade chord is 
assumed, the steady chordwise moment is 


M = EI^=— ( 103 ) 

£0 c 

where a is the coefficient of thermal expansion and AT is the difference in 
temperature between the leading and trailing edges. The corresponding static 
blade twist is 


9 


t 



(104) 
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,-6 , 0 . 


For the Sample Blade, assuming a = 15 x 10 /°F, 

G t = .225 x (500) x (1000) x(l5 x 10~ 6 )x 6.26 x 1 o‘\t 
= 1.058 x 10 _3 aT 


The temperatures of the leading and trailing edges cannot be appreciably dif- 
ferent if the chordwise camber is small. This subject will be discussed later. 

Another effect of chordwise camber is to produce a chordwise shift in the 
center of radiation pressure. The resulting pitching moment load is 


ra = p cx = 
6 n p 


g o hg o Cx p 


(105) 


where x is the distance from the center of pressure to the midchord. 
P 

dimensionless form 


„ = - P il C - X -R = 6 R 

6 p&l VcAc / P o 


( 106 ) 


Let it be assumed that the center of pressure is located at the same percent 
chord from root to tip. The blade pitch deflection for a unit dimensionless 
moment is plotted in Figure 12. Using Figure 12 and equation (106), the twist 
angle at the tip is 


e t ■ 3 - 93 ©G^ e o 


(107) 


For the Sample Blade design 


0 = 3-93 x 1000 x 6.26 x 10 






radian . 


The center of pressure shift should, therefore, be restricted to a few 
percent of chord in order to limit blade twist to an acceptable value. As in 
the case of initial chordwise curvature, the problem becomes more severe for 
longer blades. Substituting for from equation ( 1 0 1 ) gives 


4l 


= 3 - S3 G;) 


( 108 ) 


It is seen that the effect is reduced by using higher strength materials. 

Equations (102) and (108) provide practical limits to the blade radius, 


stiffness. Although it is not particularly difficult to develop concepts for 
torsional reinforcement that can be stowed in a rolled up configuration, the 
resulting increases in weight and complexity are not attractive. The analysis 
that has been presented shows that such measures are probably unnecessary for 
blades in the size range of the Sample Blade. 

A large offset in the center of pressure can also produce significant 
dynamic coupling. From the expression for the vertical component of radiation 
pressure given by case A in Table 6, the pitching moment load is 


For y = A 5 0 , ip = 90°, and 3=0. the destabilizing moment load due to 0 is 


The effect is significant, when compared with the other terms in the equation 
of motion for pitch, only when the center of pressure offset is of the order of 
several percent. Similar effects can be expected from a shift in the location 
of the center of gravity. The magnitude of the center of gravity offset can, 
however, be limited to a very small value by proper balance procedures during 
manufacture. 


which can be surpassed only by providing the blade with additional torsional 



(109) 


or, in dimensionless form, 



( 110 ) 


For the Sample Blade 


x 

rr^ = 3.76 [(0 siniji - 3 cos i|j) sin 2 y + cos 2 y] 


( 111 ) 
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The main conclusions that have been reached regarding mechanical coupling 
effects are that initial blade chordwise curvature and blade center of pressure 
shifts must be controlled. The required tolerances appear to be quite easy to 
achieve for the Sample Blade design but they become increasingly difficult for 
longer blades. 


CHORDWISE CURL 


The significance of a chordwise shift of the center of pressure relative 
to the midchord was examined in the preceding section where it was shown that 
the shift should be limited to a few percent of blade chord (for the Sample 
Blade design). An important cause of center of pressure shift is chordwise 
camber, or curl, as illustrated in the following diagram. 



The blade chord is assumed to form a circular arc of radius R c> The angle 
of the illumination with respect to a normal to the blade at its midspan is y. 
Since the radiation pressure is normal to the local blade surface, the resultant 
of the radiation pressure passes through the center of curvature and the moment 
about the blade midspan axis per unit span is 


m Q = R f 
0 c x 


( 113 ) 


where f is the horizontal component of the resultant. From the diagram 
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f =R / p sin a da 
x c 

-a /2 
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(n't) 


hk 


and, from equation (1), assuming complete reflection. 


P = P Q cos 2 (y - a) (115) 

Solution of the above three equations gives 

m e = ! p o R c(p ;) 3 sin 2 y (,,6) 

where c is the blade chord. The vertical component of force is approximately, 
for sma 1 1 a Q , 


f = p c cos 2 y 
z *o ' 


(117) 


Hence, the per unit forward shift in the center of pressure is 


x m. 

p _ 9 c tan y 

c cf &R 

z c 


( 118 ) 


The maximum pitching moment occurs when y = ^5° • for this condition x /c = c/6R . 

P c 

A reasonable design criterion for the Sample Blade desian is that x /c < 0.01. 

P 

The corresponding requirement on blade chordwi se curvature is that R^/c > 16. 67 . 

Chordwise curvature is the result of differential chordwise strain through 
the thickness of the blade, which can be produced by manufacturing procedures 
or by environmental effects. The latter, which are probably the more important, 
include thermal expansion and changes in the crystal structure of the metallic 
coating due to prolonged exposure to solar radiation. For example, solar protons 
will gradually erode a surface that is exposed to sunlight thereby producing a 
dissymmetry that will produce unsymmetri ca 1 thermal expansion, even if the 
temperature is uniform throughout the cross-section. At the present time we do 
not know the limits of the differential strain that can be produced by prolonged 
exposure in the space environment. We can, however, determine an effective upper 
limit to the moment that can be produced, as established by the yield strength 
of the material. Structural reinforcement can then be designed that will 
exhibit less than the permitted chordwise curvature while carrying the limit 
moment. 


^5 


As an example, consider the Sample Blade design which has a Mylar core with 
two external layers of pure aluminum as shown below. 



Let it be assumed that the top layer is stressed in tension to the yield 
stress and that the bottom layer is stressed in compression to the yield stress. 
The central core is assumed to be unstressed. The chordwise moment per unit 
span is 


where a is the yield stress. 

y 

The moment produced by differential expansion of the surface layers is 
resisted by the bending stiffness of the core material and by the chordwise 


where x is measured from the mid-chord. It is easily shown that neither effect, 
is sufficient to reduce the chordwise curvature to an acceptable level and that 
reinforcing members are required. 

Figure 13 shows a design concept for a chordwise reinforcing member, or 
batten, that was developed at Astro Research Corporation as part of the effort 
described in this report. Each batten consists of a thin film of Kapton on 
one side of which a carbon coating has been deposited. Chemical interaction 
between the Kapton and the carbon produces a large prestress which, upon release 
of one dge of a strip of the coated film causes it to roll into a tight scroll 


m = h h, o 
c f y 


(119) 


component of centrifugal force which produces a chordwise tension stress (Ref. 9), 



( 120 ) 
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of approximately 1/8 inch diameter. A volatile material which diffuses rapidly 
in vacuum provides a temporary bond on one edge of the strip so that it may 
conveniently be rolled in a flattened condition into the stowage roll. 

Structural test of a sample scroll showed that it had a bending stiffness 

(El)^ equal to 6.39 lb-in 2 and that its maximum bending moment capacity, as 

limited by local crippling, was 0.265 lb-in. The mass per unit length of the 

-A 

scroll was 1.46 x 10 lb/in. 


It is important to attach the batten to the blade in such a way as to per- 
mit differential expansion in the chordwise direction between the blade and the 
batten. This can be done by slotting the batten scroll near the edge that is 
attached to the blade. If the scroll is slotted, the total bending stiffness 
is the sum of the stiffnesses of the blade and of the batten. Since the latter 
is much larger than the former, the chordwise bending moment developed per unit 
span is 


m 


c 


(El), 


R c*b 


( 121 ) 


where is the distance between battens. Thus, from equations (119) and (121) 


(El). 

^b R h h, a 
c f y 


( 122 ) 


For the Sample Blade, using the batten design concept that has been described 

and the design criteria that have been suggested, the parameters in equation (122) 

are 


(El) b = 6.39 lb-in 2 

R = 16.67 c = 2000 in 
c 

h = 2.5 x 10 ^ in 

h = . 15 x 10 ^ cm = .0591 x 10 ^ in 
0^ = 5000 psi (pure aluminum). 


The required batten spacing is 


*b = 


6.39 


V 


2000 x 2.5 x 10 x .0531 x 10 x 5000 


= i )32 in = 36.0 ft. 


The bending moment carried by each batten is M = (El) . /R = .0032 in-lb, 

c 13 c _/j 

which is well below its capacity. The weight of each batten is 120 x 1.46 x 10 = 

.01753 lb. The required number of battens is 10,000/36 = 278. The total weight 
of all battens is .01753 x 278 = A. 88 lb, which is 2.5% of the blade weight. 


Since the batten spacing is large compared to the chord, the manner in 
which blade deformations are restrained by the battens should be investigated. 
This may be done by examining the following differential equation which governs 
the deflection of the blade normal to its surface. 
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where 


h/2 



-h/2 


dz 


( 12 *)) 


is the bending stiffness. Since the curvature in the spanwise direction is 
small, terms © and © may be neglected. Term (j) is very small, because a x 
is very small (see equation (120) and may be neglected. 


The boundary conditions are that the vertical displacement is restrained 
by battens at r = r^ and at r = r^ + and that the leading and trailing edges are 

free. The normal pressure, p n , is very small, but differential expansion of the 
blade produces an initial curvature whose effect may be represented by placing 
bending loads equal to m^ on the leading and trailing edges. For simplicity, 
it will be assumed that the battens are undeflected and that the deflection of 
the blade varies sinusoidally in the spanwise direction. Thus 
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w(x) 

m 
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where Ar = r - With this assumption, and omitting the terms that have been 

declared negligible, the differential equation is 


8 4 w 
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ax 4 


+ (—■'f a h w = 0 

VV r m 


( 126 ) 


At the leading and trailing edges, x = ±^-, the boundary conditions are 

8 2 w m_ 
m _ c_ 
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3x 3 


(127) 


= 0 


where, again, curvature in the spanwise direction has been neglected. Equations 
(126) and (127) constitute a two-point boundary value problem that is analogous 
to the problem of a uniform beam with end moments on an elastic foundation. The 
solution to the problem is 
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where 


tt 2 o h\l/4 
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(129) 


etc 

For Heliogyro blades, is a large number which results in simplification. For 
the Sample Blade 


D = 2.75 x 10“ 6 lb-in 

a r = 750 psi (at midspan) 

- 

h = 2. 5 x 1 0 in 

H b = 432 in 

c = 120 in 
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so that 


a = 1.38 in 


and 


f- = 82.9 


Thus, in equation (128) ctnh = 1.0 and cosh ax/sinh li- 
near the edges. The resulting simplified expression for the 


« 1 except very 
deflection is 
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which shows that the chordwise deflection shape is a wave. The chordwise wave- 
length X = 2-nVa is equal to ^.56 inches for the Sample Blade. The wave height 
i s 


w. = 


cos 


ac 
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+ s in 


ac 
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(131) 


Since y- is a large number, the wave height can change a large amount for a 
small change in a. The maximum value is 


hm Dcr 


(132) 


For the Sample Blade, m^ = h h^. 0 ^ = 7-39 x 10 ^ in-lb/in, and 


hm 


7.39 x 10 
■?> 


-6 


= 1 .1*1 in. 


2.75 x 10 0 x (1.38)' 


It remains to examine whether a wave of this height produces a significant shift 
in the center of pressure. It may be noted in this connection, that the wave 
height is directly proportional to batten spacing, as may be seen from equations 
(129) and (132). Thus the wave height can be reduced by reducing the batten 
spacing. 
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An analysis of the radiation pressure acting on a wavy surface whose 
deflection is 

ac 


w = w. 


s i n 
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3w 


sinh 


cosh ax 


( 133 ) 


produces, under the assumption that — « 1, the result that the pitching moment 

oX 

per unit span about the mid chord is 


= cp o w h s i n 2y [cos - sin |£j 


(13A) 


Note that the magnitude of the moment fluctuates between ±cp Q w^ sin 2y for small 
variations in the physical parameters because ac/2 is a large number. This 
characteristic may, in itself, be relied upon to produce an average pitching 
moment load over a relatively long blade span that is much smaller than its 
maximum value. The maximum value of m is 
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(135) 


The per unit shift in the center of pressure is, from equations (117) and (135) 


x m~ w, 

= __ = 2 — tan y 
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(136) 


Thus, for the Sample Blade and the condition of maximum moment (y = ^5°), 

.0333 


2/?x '- 1 " 


120 


Considering that a very severe differential strain condition (yield stress on 
top and bottom surfaces) has been assumed and that the effects of the variations 
in physical parameters have been ignored, this result shows that the center of 
pressure shift under practical conditions is probably negligible. A more im- 
portant concern may be the presence of deep spanwise ridges in the blade which 
will produce chordwise shrinkage and which will tend to invalidate the assumption 
of specular reflection which is made in performance calculations. A shorter 
distance between battens will tend to alleviate the problem, but a better solu- 
tion is to establish a less conservative estimate of differential expansion by 
further research. 
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DYNAMIC RESPONSE DURING MANEUVERS 


A major concern in the design of a flexible vehicle is the extent to which 
structural flexibility will interfere with the response of the vehicle to 
maneuver commands. This is particularly true of the Heliogyro because of the 
extreme flexibility of its blades. Accordingly, it was deemed important to 
perform an analysis of the dynamic response during maneuvers under the present 
study. 

The primary means for exerting forces on the vehicle in the plane of rota- 
tion is cyclic blade pitch, as shown in Figure l^a. A simple application is 
shown in Figure 14b, where cyclic pitch is used to change the orbital altitude 
of a vehicle in a polar orbit. Note that the frequency of the cyclic pitch input 
should be equal to fi - \p, where \p is the orbital frequency, in order to align 
the resultant control force with the orbital velocity at all times. 

The theory for calculating the longitudinal and lateral forces due to cyclic 
pitch is relatively simple if it is assumed that the blade is rigid. Consider, 
for example, the case of axial illumination and small pitch angle. Then, the 
total force on the blade in the chordwise direction is, from Table 6, case A, 


F = p Rc0 
x K o 


037 ) 


The force components in the longitudinal and lateral directions (x and y 
in Figure 11) are 



( 138 ) 


F, = -F cos fit 
lat x 


Let the blade pitch be cyclicly excited at frequency fi + u) q , such that 


0=0 sin (fi + 0 ) ) t 
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Then, combining the above equations and resolving the products of trigonometric 
functions into sum and difference terms, 


F, = ~ p Rc0 [cos to t - cos (2ft + (jj ) t] 
long 2 o c o o 

F, = — »■ p Rc0 [sin at t + sin (2ft + a) )t] 
lat 2 o c o o 


0*tO) 


The high frequency terms have little significance and they will, in any 
event, be cancelled out when the effects of three or more blades are added. The 
low frequency terms are the desired control forces. 

The elastic response to the sudden application of cyclic pitch was analyzed 
during the present study by means of a digital computer simulation. The pitch 
input was assumed to be a sine wave which is modulated by an envelope function 
that rises to its full amplitude in one-half revolution, see Figure 15a. Three 
different cases were analyzed corresponding to different values of the frequency 
of the sine wave. Blade pitch was assumed to produce a component of force in 
the inplane direction according to equation (137). which is the appropriate 
expression for the case of axial illumination and small pitch motions. 

The simulation included a single blade, with properties identical to those 
of Sample Blade design. The center body was assumed to be massive compared to 
the blades so that the translational motions at the blade root could be rigidly 
restrained in the simulation. It is evident from the symmetry of the applied 
load, see Figure l^a, that the inplane bending slope, t , , should also be rigidly 
restrained at the blade root. 

Mechanical coupling between blade motions was ignored in view of the demon- 
stration (see section on MECHANICAL COUPLING DUE TO BLADE DEFLECTIONS) that these 
effects can be made negligible for the Sample Blade by careful construction. As 
a result vertical deflections could be ignored entirely since they were not 
coupled, either directly or indirectly, to the source of dynamic excitation. 

Dampers were included at the blade root for torsional and inplane bending 
motions. The properties of the dampers that were used and the damping that 
they provided for the lowest four modes of each type are given in Tables A and 5* 


53 


The 

ana 1 y 1 

tical 

mode 1 i 

spanwi se 

pos i t i 

i ons : 


STA 


1 

2 

y = 

r/R 

.005 

.010 

Note the 

concentrati 

ion of s 

adequate 

s i mu 1 ; 

ation 

for the 


ncluded ten blade stations 

3 4 5 6 

.020 .040 .080 .150 

tat ions near the blade root 
rapid decay of chordwise b 


ocated at the following 

7 8 9 10 

.30 .50 .70 .90 

in order to provide 
inding moment. 


The results of the analysis are presented in Figures 15a thru 152, for the 
case in which the frequency of the cyclic pitch excitation is equal to the 
rotational speed. The response of various quantitites are plotted for the first 
twenty revolutions following the initiation of the maneuver. The plots were 
produced by a computer, which accounts for the somewhat unconventional scaling 
practices that are evident. 


The blade pitch response is shown in Figures 15b, c and d, where it is seen 
that the pitch motion achieves its steady state value (within ±10%) in about five 
revolutions. 


The inplane blade deflections at 50% and 90% span are shown in Figures 15 e 

and 1 5 f • The presence and the decay of the low frequency (.105/rev) inplane mode 

are evident in these plots. Note also that the amplitude of motion at the tip is 

comparable in magnitude to the static vertical tip displacement (.694 R0 ) when 

o 

the magnitude of the pitch input is one radian. This is the order of magnitude 
of dynamic chordwise motion that was assumed earlier in estimating mechanical 
coupl i ng effects. 

The inplane and axial components of force that exist at the blade root are 
shown in Figures 15g and 1 5h . The axial component was obtained by summing the 
axial component of Coriolis force, f in equation (89), along the blade span. 

The inplane component is the sum of the inplane bending shear and the inplane 
component of centrifugal force at the blade root. It is seen that the axial 
component is, surprisingly, about twice as large as the inplane component. 

The forces at the blade root are resolved into components parallel to non- 
rotating axes in Figures 1 5 « and 1 5 j - Note the presence of the steady component 
of longitudinal force and the existence of the 2/rev oscillations as predicted 


5 ^ 


by rigid blade theory. It is seen that the low frequency component of the 
longitudinal force settles down to the neighborhood of its final value within 
about five revolutions. The gradual decay of the high frequency transients can 
also be detected. 

The inplane bending moment near the blade root is plotted in Figures 15k and 
15&, where it is seen that the bending moment decays to 10% of its value at the 
root in only 2% of the span. The peak magnitude of the inplane moment at the 
root is 


M = .006 p R 2 c 9 
l, o c 

= 13-5 6 in-lb 
c 

The corresponding displacement of the tension axis, as a fraction of blade chord, 
i s 

x M '3-5 e c 
f - TS ” 30TT20 ’ •° 0375 e c 

Thus, even a one radian pitch amplitude produces a shift in the tension axis 
that is negligibly small. 

The longitudinal and lateral force resultants produced by the simulation 
study are compared with the results of rigid blade theory in Figure 16 for three 
values of the cyclic pitch frequency, w. The 2/rev components of force were 
filtered from the results of the simulator study. It is seen that, in all cases, 
control is established within approximately one revolution and that the steady 
state is achieved (within ±10%) in about five revolutions. It will be noted 
that the steady state forces for the flexible blade exceed those for the rigid 
blade when the pitch frequency is greater than one/rev, and are less than those 
for the rigid blade when the pitch frequency is less than one/rev. This result 
is explained by Figure 5. 
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EXPERIMENTAL PROGRAM 


An experimental program, which was designed to demonstrate the deployment 
concept and the control of blade pitch, was conducted at Astro Research Corpor- 
ation as part of the present study. It was realized in advance that it would 
be impossible to simulate the proper values of the critical dynamic parameters 
(K, 3 q ) in an Earth environment. At the same time it was hoped that the exper- 
iment might provide a verification of the theory of dynamic response. As 
matters turned out the experiment was, in fact, a critical test of the theory 
for dynamic coupling effects due to blade deflections. 

The experimental apparatus consisted of an 18 Ft. diameter spin chamber in 
which a model of the Heliogyro blade was deployed from a spool. The spin chamber 
was not evacuated, but since it was enclosed and rotated at the same velocity 
as the experimental model, the velocity of the enclosed air relative to the 
blade was substantially zero and the aerodynamic effects (virtual mass and damp- 
ing) could be neglected. 

A diagram of the test set up and the parameters that define the experiment 
are shown in Figure 17* It will be noted that the deployment spool is offset a 
considerable distance from the axis of rotation. Note also that the pitch axis 
of the blade is tangent to the surface of the spool and that it is inclined 
downward at an angle that is somewhat less than the root coning angle of the 
fully deployed blade. A number of experiments were conducted including unwind- 
ing and rewinding the blade, steady cyclic pitch oscillations, and abrupt 
changes in collective pitch. A motion picture was produced that demonstrates 
successful accomplishment of these maneuvers. 

The experimental blade consisted of a .001 inch thick strip of Kapton about 
2 inches wide and 77 inches long. Photographs of the deployed blade made from 
motion picture frames at two different pitch angles are shown in Figure 18. The 
inplane deformation of the blade is clearly visible in the lower photograph. 

The blade was painted in order to enhance its visibility and to provide lateral 
markings for the observation of pitching motions. The paint, which was applied 
to one side only, also produced chordwise camber. 
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The frames of a motion picture showing the response of the blade to an 
abrupt change in collective pitch from zero degrees to twenty-five degrees were 
analyzed to produce a time history of the pitching motion at the 50% span and 
at the blade tip. The results of a digital computer simulation are compared 
with the experimental measurements in Figure 19. The agreement between analysis 
and experiment is seen to be only fair. I t wi 1 1 be shown that the most likely 
reasons for the observed discrepancies are nonlinear effects, not accounted for 
in the analysis, and the extreme sensitivity of the response to the angle of 
depression of the pitch axis. 


The most important dimensionless modeling parameters for comparing the 
dynamic response of different blades are the inplane bending stiffness parameter, 
K, and the steady coning angle at the blade root, 8 q . The values of these para- 
meters for the Sample Blade and the experimental blade are compared below. 


Sample 

Exper i menta 1 

Experimental Blade 

Blade 

Blade 

Sample Blade 

.833 x 10"* 

. 45** 

5440 

6.26 x 10 * 

.1967 

314 


The conclusion stated earlier, that dynamic coupling due to blade deflec- 
tions is unimportant for the He 1 iogyro, obviously does not apply to the experi- 
mental blade. The coupling terms in question are written in dimensionless 
form in equations (93) , (94) and (95). If small dynamic motions are assumed, 
the important terms are those that contain the static coning angle = 3w s /3y. 
They are 


where S 1 


ae s 92w s 

“ay = * 


f 


x 


E(6 S ') 

2a 

o 

f = 


^s 3T + 


- e + 

p s 3c 

(141) 

2o o 9 y 

i is. 

s 9T 


(142) 

1 9 2 

2 3 r 

r (K$ s 'e) 

(1**3) 
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The analysis was performed including the above terms. The boundary condi- 
tion at the blade root accounted for the angle between the blade twist axis and 
the applied pitch axis. Thus the inplane blade slope at the root is 

c o = l7 = _0 c sin (b o " y 


The sensitivity of the analysis to the inclination of the pitch axis, Bp, 
is indicated in Figure 20 where it is seen that a change in B of one degree 
produces a 20 degree change in the peak pitch angle at the tip. 


Since the peak transient pitch angle at the tip is actually rather large 
(72°) the error produced by the small angle assumption in the analysis cannot be 
neglected. Furthermore, the peak magnitude of the inplane bending moment at 
the blade root is large enough to produce buckling at the edges of the blade, 
thereby altering the inplane bending stiffness of the experimental model signi- 
ficantly. 

The effects of finite elastic torsional stiffness and of the added mass 
due to the air were investigated and were found to be unimportant. The most 
probable reasons for the differences between the analytical and experimental 
results are the effects of large dynamic motions, the effects of large initial 
deflections, and the sensitivity of the response to the inclination of the 
pi tch axis. 
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CONCLUSIONS 


A number of fundamental topics concerning the structural dynamics of the 
Heliogyro have been examined. These topics include 

— Uncoupled vibration modes 
— Means for providing damping 
— Coupling between the blade motions 
— Transient response to maneuver commands. 

On the basis of the results of the study it may be concluded that: 

1. Resonances between the harmonics of rotor speed and the blade vibration 
modes are easily avoided. 

2. Damping can be provided for all degrees of freedom. The damping in 
pitch is sufficient to produce convergence to within 10% of steady 
state in about five revolutions. The damping that can be provided 
(with very little added weight) for inplane and vertical vibrations is 
about an order of magnitude smaller but it is probably adequate. Damp- 
ing for the rigid body nutation mode requires active control feedback 
or an auxilliary passive damping device connected to the center body. 

3. The major source of coupling between the blade motions is solar radia- 
tion pressure which produces inplane and vertical forces due to pitch- 
ing motion. This source is utilized in orientation and control of the 
vehicle. 

4. Pitching moments due to vertical and inplane blade deflections can be 
held to an unimportant level for blades that are the size of the Sample 
Blade (10,000 ft radius) by careful design and manufacture. The major 
sources of unwanted pitching moments are a chordwise shift in the loca- 
tion of the center of pressure due to chordwise camber, and an initial 
chordwise radius of curvature in the as-manufactured blade. Suppression 
of chordwise camber requires chordwise bending members (battens) at 
frequent intervals. The design of a satisfactory batten which increases 
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the weight of the Sample Blade by 2.5 % has been demonstrated. The 
control of initial chordwise curvature to the level required by the 
Sample Blade does not appear to be difficult. 

5. Undesirable dynamic coupling effects (pitching moments) increase in 
severity with increases in blade radius. For example, the steady 
pitch angle at the blade tip due to center of pressure shift, x^, 
is, from equation (108), 


6 


t 





and the steady pitch angle at the tip due to initial chordwise radius 
of curvature, R^, is, from equation (102), 



It is seen that in both cases the pitch angle increases in proportion 
to the square of the span, R, and inversely as the thickness, h. The 
remedies, in addition to suppressing the primary sources, x^ and 
1/R c , are to increase the blade chord, c, and, more importantly, to 
increase the blade centrifugal stress, a Q . Thus, the development of 
higher strength materials will enhance the feasibility of designing 
longer blades. Another remedy, which has the undesirable effects of 
increasing the weight and complexity of the blade, is to provide 
torsional stiffness by means of structural reinforcement. 


6. The effect of blade flexibility on the ability to provide lateral 

control forces by means of cyclic pitch is slight. Effective response 
to an abrupt command is achieved (within ±10%) in about five revolu- 
tions. The effect of rotor size on control effectiveness is small. 


In summary, it is concluded that no structural dynamic reasons to doubt 
the feasibility of the Heliogyro in sizes up to 10,000 ft radius have been dis- 
covered. On the other hand, the development of very much larger sizes will 
depend, in part, on the ability to minimize undesirable dynamic effects. 
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APPENDIX 


THE TWISTING MOMENT IN A SLIGHTLY CURVED BEAM 

Consider a beam whose axis is slightly deflected from the r-axis, as shown 
in the following sketch. 


z 



The angles and £ , which define the orientation of the pitch axis at the root 
end, are small. Thus, to first order in the angles, the moment about the pitch 
axis is 


M = M + 
a r 


£ M + 
o x 


e m 

o z 


(A-l) 


The moments about the x, r and z axes due to running loads applied to the 
centroidal axis of the beam are, from elementary mechanics, 
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■■■/ 
o 

"*■ / 


z + f r)dr 
z 


(f z - f x)dr 
x z 


(-f r + f x) dr 
x r 


Upon substitution into equation (A - 1 ) 

.R 




[f z - f x + C (-f z + f r) + 3 (-f r + f x)]dr 

x z o r z o x r 


(A— 2 ) 
(A- 3) 
(A-A) 


( A- 5 ) 


Assuming that x, z, f , f , and C and 3 are of first order and that f and r 

are of zero order, it is seen that all of the terms in equation (A-5) are of 

second order. Equation (A-5) is inconvenient for use in an analysis where the 

other effects are expressed by partial differential equations because of the 

presence of 3 and £ . It can, however, be converted into a more amenable form 
o o 

by applying the rules for integration by parts. First, however, it is convenient 
to rewrite equation (A-5) in terms of the offset distances 


z = z - 3 r 
o 


x = x - C r 
o 


(A-6) 


from a line tangent to the pitch axis at the root. Thus, 

r R _ 

M a = / [f x z ' f z x + f r (x3 o " zC o )]dr 


(A-7) 


The formula for integration by parts is 


f b 1* 

I udv = uv - / vdu 


(A-8) 


The formula will be applied in the following manner to the terms in equation (A-7) 
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Te rm 

u 

dv 

du 

V 

f 7 

X 

z 

f dr 

X 

(3 

- 3 ) dr 
o 

f R 

- / f dr = -F 

J X X 

r 

f X 

z 

X 

f dr 
z 

(C 

- ?Q )dr 

- f R f dr = -F 

l 2 

f (xB - z c) 
r o o 

o o 

f dr 
r 

O 

- BC 0 )dr 

- j* f r dr - -T 
r 


The quantities F x> and T that appear above are respectively the total inplane 
force, the total vertical force, and the total axial tension at a distance r 
along the blade. The result of the integration by parts is 


M 

a 


zF + xF 
x z 


- T(x6 o -z ? o ) 


R 

+ 

o 



[F x (8 ‘ S o ) 


F z <?-?„> + T(t;S o -6c o )]dr (A-S) 


The boundary terms vanish because z, x, x, and z are zero at r = 0, and 

T are zero at r = R. The inplane force, F , and the vertical force F 

x z 

carried by centrifugal tension and partly carried by beam shear. Thus 


F , F and 
x z 

are partly 


and 


F = T£ 



F 

z 


= T6 


8 r 


(A- 10) 
(A- 1 1 ) 


where M and M are respectively the inplane and vertical bending moments. Upon 
s P 

substitution of equations (A-10) and (A- 1 1 ) into equation (A-9) , it is seen that 
all of the terms proportional to T cancel each other so that 


M 


-/v 


3M. 


3M 


^(B-6 0 ) + ^(c-g)dr 


(A-12) 


Integration of this equation by parts produces the result 

-R 


"a ’ /<"; I? - " 6 £ )d ' 


(A- 13) 
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Since the origin can be shifted without changing the form of the result, it is 
seen that the running pitching moment load is 


m e = \ 97 


M *1 
0 3r 


(A- 14) 


The term indicated in equation (A— 1 ^ ) can be added directly to the dif- 
ferential equation for pitch. 
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TABLE 1 


SOLAR SAIL DESIGN MORPHOLOGY 


SHAPE 

MEANS OF RIGIDIZING 

MEANS OF ORIENTING 

CIRCULAR 

ELASTICITY 

AREA SHADING 

'J/ 



RECTANGULAR 

ELECTROSTATIC FORCE 

CENTER OF GRAVITY MOVEMENT 

'i' 



VERY SLENDER 

MAGNETIC FORCE 

REFLECTING SURFACE ROTATION 


CENTRIFUGAL FORCE 



; 
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TABLE 2 


PROPERTIES OF SAMPLE BLADE 

DESIGN 

Blade Chord (uniform) 

10 ft. 

Blade Radius 

10,000 ft. 

Blade Material 

Chordwise Bending Moment of Inertia 
Surface Density 

Modulus of Elasticity (Composite) 

1 /A mi 1 My 1 ar wi th 
1500 A° aluminum 
coating each side. 

36 i n 1 * 

1.3^ x 10 ps i 

1.0 x 10 6 psi 

Rotational Speed, ft 

.0316 rad/sec 

Rotational Period 

3.31 minutes 

Blade Root Stress-- 

1000 psi 

Blade Root Tension 

30 lbs. 

B 1 ade Weight 

1A lbs. 

Total Thrust for Normal Solar Illumination 

.0188 lb. 

Sail Lightness Number, 
Blade Root Coning Angle 

0.16 

6.26 x 10 4 rad. 

* Based on lA mil thickness. 
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TABLE 3 


UNCOUPLED BLADE FREQUENCIES 


u = w/ft 

n 

TWIST 

VERTICAL 

INPLANE 

i 

YiT = 1.414 

1 .0 

0 

3 

yy = 2.646 

yr" = 2.449 

VT = 2.236 

5 

y y = 4.000 

yry = 3.873 

yy = 3.742 

7 

yiy = 5.385 

yyr = 5.291 

yy = 5.099 
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TABLE 4 


FREQUENCIES AND DAMPING FOR BLADE 
PITCH MODES WITH A CONTROL DAMPER. 


K d 

pft 2 !R 


2 0 -—= — 
’ k , n 


Mode 

Number 

Frequency 

Cycles/Rev. 

Damping 

9 

Revolutions to 

1 /2 Amp 1 i tude 

1 

1.368 

.0845 

1.91 

2 

2.618 

. 1024 

.825 

3 

4.022 

.0863 

.636 

4 

5.455 

.0680 

.596 


TABLE 5 

FREQUENCIES AND DAMPING FOR INPLANE 
MOTIONS OF THE SAMPLE BLADE WITH A DAMPER 
WHOSE PROPERTIES ARE GIVEN IN EQUATION 65. 


Mode 

Number 

Frequency 

Cycles/Rev. 

Damp i ng 

g 

Revolutions to 
1 /2 Ampl i tude 

1 

.1053 

.1167 

18.0 

2 

2.276 

.00236 

41.0 

3 

3-834 

.00196 

29.5 

4 

5.342 

.00159 

26.1 
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TABLE 6 


FORMULAS FOR VERTICAL AND INPLANE LOADING DUE TO SOLAR ILLUMINATION 


ASSUMPTION 

f 

rz 

f 

rx 

CASE 

6 

0 

Y 

1. 

0 

0 

any 

p c cos 2 y 
o ' 

0 

2. 

0 

« 1 

any 

P q c[cos 2 y - 6 sin 2y cos 

0 

3. 

any 

any 

0 

p c cos 2 6 cos 3 0 
o 

p c cos 2 0 cos 2 0 sin 0 
o 

A. 

« 1 

« 1 

any 

P q c[(0 sin \p. - 0 cos ip) s i n 2y + cos 2 y] 

p c 0 cos 2 y 
o ' 

5. 

any 

« 1 

any 

P q c cos ©[(sin 0 sin sin y + cos 0 cos y) 2 

-0 sin y (^- sin 20 sin sin y 

+2 cos 2 0 cos \p cos y)l 

f rz tan 0 


> 




Fi gure 


Solar Radiation Pressure on a Reflecting Surface 
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1 , 000,000 


Spec i f i c 
Impu 1 se 
(sec. ) 


100,000 


10,000 


1,000 


100 



10 


10 


100 

Mission Duration 


1,000 

(Days) 


10,000 


FIG. 2. SPECIFIC IMPULSE VS MISSION DURATION FOR VARIOUS PROPULSION 
SYSTEMS OPERATING IN THE VICINITY OF THE EARTH'S ORBIT 
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Figure 3. Sketch of Experimental Two Blader 
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Figure b. Rotor Blade Coordinate Notation 


7b 


Figure 5. Blade Tip Pitch Response due to 
Sinusoidal Excitation at the Blade Root 



u _ Pitch Frequency 
Rotational Speed 
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BLADE 




Figure 6. Idealized Model for Blade Pitch 


t 



Mg = Ay; K 0 = Ay; K 


Pitch Angle at Blade Root 



Pitch Angle at Blade Tip 



Figure 7* Blade Torsional Response to Ramp Collective Input 
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Figure 8. Results for Control Damper Optimization Study 



(a) 2nd Collective Inplane Mode 



(b) 2nd Cyclic Vertical Mode 


Figure 9- Blade Modes that Result in 
Rotations of the Center Body 
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8 



A-A' 


Dr i vs 
Shaft 



Figure 10. Blade Retention System 


z 



Figure 11. Coordinate Geometry for Rotating Blade 
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Figure 12. Pitch Distribution Due to Static Applied Load 
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73 |-| 


Temporary bond 


Permanent bond 


Stowage rol 1 


Figure 13- Deployable Transverse Batten Concept 



(a) Lateral Force Produced by Cyclic Pitch 

0 = 0 s i n ftt 
c 



(b) Use of Cyclic Pitch to Change Orbital Altitude in a Polar Orbit 
Figure 1**. Lateral Vehicle Control 
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(b) Pitch Angle at Blade Root 
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00 

Revs . 
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2.00-1 


(c) Pitch Angle at 50% Span 
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00 


Revs . 
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(h) Axial Force at Blade Root 
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(X10 l ) 
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Figure 15. Response to Cyclic Pitch Input 
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OQ 
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(X10"®) (k) Inplane Bending Moment 



(X10 ») 

00 Revs . 


( X 1 0 ) (£) Inplane Bending Moment 
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Longitudinal Force 


(a) w = ft 


Lateral Force 





a) = 0.9Q 




Figure 16. Force Resultants for Cyclic Pitch Input 

Flexible Blade 
Rigid Blade 
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ft 



Blade Properties 


Thickness: .001" (Kapton) 

Chord, C = 1.938" 

Weight/unit area = 1.08 x 10 ^ lb/in 2 (Kapton plus paint) 
Modulus of elasticity, E = 500,000 ps i 

Operating Parameters 

Spin rate, ft = 2 tt rad/sec 

Stress at Centerline, a = ^3.6 ps i 

o 

Coning angle at tip, 0 t = 6.33° 

Inplane Stiffness Parameter, K = jy (a~ )(^") = 


Figure 17- Heliogyro Blade Experiment 
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Figure 18. Photographs of Rotating Blade in Test Chamber 
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(a) Pitch Angle at 50 % Span 



(b) Pi tch Angle at Tip 
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6=7° 

P 



3 P = 8° 



3 P = 9 ° 



Figure 20. Effect of Pitch Axis Inclination on 
the Dynamic Pitch Response at the 903 Span 


NASA-Langley, 1971 32 
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